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We extend the stochastic energetics to a relativistic system. The thermodynamic 
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relation between the relativistic equipartition relation and the rate of heat transfer 
is discussed in the relativistic case together with the nature of the noise term. 



Q^. I. INTRODUCTION 

D 

Physical interpretation and theoretical interpretation of Brownian motion, a well-known 

(N : 

>• . process named after the Scottish botanist Robert Brown in 1827, were first established by 

in 

00 ' Einstein and Smoluchowski, then developed by Langevin and many others. From these 

m ■ 

■ studies, it was made clear that Brownian motion is nothing but the manifestation of the 
Q , presence of invisible microscopic bodies such as atoms and molecules, and their thermal 
motion. Thus, Brownian motion has naturally been discussed in relation with kinetic theories 
and thermodynamics. Nowadays concepts and methods of stochastic processes introduced 
in the formulation of Brownian motion are widely applied to various fields of science such 
as biophysics and economy. Mathematical foundation of these methods has also been well 
studied. 

After thermodynamics was established in the nineteenth century, many efforts have been 
made to understand the thermodynamic principles from the microscopic point of view. The 
achievements of these efforts are summarized today as statistical mechanics. There, the time 
average of microscopic behaviors is replaced by an average in terms of a suitably chosen 
statistical ensemble of microscopic events and the temporal behavior of each particle is not 
questioned. Therefore, the basic question to the foundation of statistical mechanics is how 
this substitution of averages is justified from the microscopic deterministic point of view. 
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Of course such a question has been studied deeply in many ways since the time of Maxwell, 
Boltzmann and Gibbs such as if-theorem and ergodic theory. Mathematically rigorous or 
not, for physicists the overwhelming success of statistical mechanics itself to describe the 
nature for thermally equilibrated systems could be considered as the proof of validity of the 
hypothesis. 

On the other hand, recently, it is more and more becoming important to understand 
the bulk properties of a dynamical system where the thermodynamic equilibrium is not 
necessarily attained. For example, in the field of relativistic heavy ion collisions (RIG), we 
would like to determine the thermodynamic properties of quark and gluon plasma (QGP) 
from the experimental data. RIG is somehow the unique way in laboratories to extract the 
thermodynamic properties of Quantum Ghromodynamics (QGD) in the deconfined phase. 
However, of course, real systems in RIG situation are not in thermodynamic equilibrium, 
since the system size and the reaction time scale are finite. Although many successful 
models, such as relativistic hydrodynamics have been employed to describe several aspects 
of the observed data, we feel still an existence of "missing link" between the thermodynamic 
properties and microscopic dynamics. For example, several questions like, "How much can we 
be sure that the equation of state used to fit the experimental fiow data with hydrodynamics 
should be the same as that of obtained from the lattice QGD calculations?" , "Is the extended 
thermodynamics valid and unique?", or "How to deal with the relativistic covariance and 
thermodynamic limit?" are yet difficult to be answered. In this sense, it may still worthwhile 
to ask the repeated question for the basic foundation of statistical mechanics from the view 
point of these new physics necessities of nowadays. 

The stochastic energetics may serve a useful tool to understand the questions above 
In this approach, the dynamics of each event of a Brownian motion is directly refiected in 
the change of thermodynamic quantities, such as heat, energy and work. That is, these 
thermodynamic quantities are also dealt as stochastic variables, refiecting in the name of 
this approach, stochastic energetics. Therefore, the application of this approach seems to be 
useful for RIG studies, where event by event fiuctuations are considered to carry important 
information on QGD [2]. 

So far, the stochastic energetics has been applied to non-relativistic phenomena. However 
it seems very attractive to apply this approach to investigate the thermodynamic properties 
for relativistic systems. This is because it is possible to incorporate the effect of special 



3 



relativity (such as the Lorentz transform) into the model in a straight-forward manner, per- 
mitting to discuss directly the relation between thermodynamics and the theory of relativity. 
Such an analysis is difficult in statistical physics where temporal averages are replaced by en- 
semble averages from the beginning. The purpose of this paper is to show that the stochastic 
energetics is still applicable to Brownian motion of a relativistic particle. It is important 
to notice that this is not a trivial problem, because equations of motion for stochastic vari- 
ables involves a finite time interval so that the space and time entangles when the Lorentz 
transformation is involved. Therefore, we have to pay a special attention for the treatment 
of noise when we discuss a relativistic stochastic motion in a general reference frame. 

The present paper is organized as follows: In Sec. 2, we introduce a model of relativistic 
Brownian motion. Relativistic generalization of Brownian motion has been discussed by 
many authors {s-lSj, but the definitive answer to this problem has not yet established, 
or simply does not exist in a unique manner. As mentioned above, this is because of the 
difficulty to introduce a manifestly Lorentz covariant noise in an unambiguous way, as well as 
its physical model. Here, instead of formulating a covariant equation of motion for stochastic 
variables in a general case, we discuss a less ambitious problem. We start from a Langevin 
equation for relativistic Brownian motion defined in the rest frame of a heat bath. Then 
we determine how the noise term should transform under an arbitrary Lorentz boost of the 
system in order to keep the internal consistency. We show that in such a situation, we can 
introduce a consistent noise term. In Sec. 3, we apply the idea of the stochastic energetics 
to our model of relativistic Brownian motion and derive the first and the second law of 
thermodynamics. We also discuss the relation between the rate of heat transfer and the 
relativistic equipartition relation. Summary and concluding remarks are given in Sec. 4. 

In this paper, we used the natural unit, h = c = ks = 1- 



II. MODEL OF RELATIVISTIC BROWNIAN MOTION 

As mentioned in the introduction, the formulation of the relativistic Langevin equation 
contains still an open question due to the transformation property of the thermal noise. To 
avoid this question, we start with a well-defined model of relativistic Brownian motion in 
the rest frame of a heat bath. Then we study the transformation of the system under a 
Lorentz boost and check the consistency between our model and special relativity. 
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A. Relativistic Brownian motion in the rest frame of heat bath 

We consider the Brownian motion of a relativistic particle with mass m in the 3+1 
dimension. In the rest frame of a heat bath, we consider the following Langevin equation, 

^ = (la) 

dt* pO* ^ ^ 

^ = -K/*)p* + Vm^^it*), (lb) 



where p'^* = a/ {p*y + m^, that is, the particle always kept on the mass-shell. Throughout 
this paper, we use the symbol * to indicate the value of a variable defined in the rest frame of 
the heat bath. The parameters v{jp*) and D{p^*) characterize, respectively, the relaxation 
of the momentum and the strength of the noise. We assume that they are Lorentz scalar 
functions, depending only on the particle energy in this frame. N is a Gaussian white noise 
three- vector 



and has the following correlation properties. 
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(2) 



(N(r))RF = 0, 

{N\t*)N^it*'))KF = 5,,5it*-t*'). 



(3) 
(4) 



The symbol (X)r,^ denotes the stochastic average of X in the rest frame of the heat bath 

nn 

(we refer to as RF). The similar Langevin equations have already discussed in Refs. |10|-|15|. 
Note that the Langevin equation can be obtained even from a binary collision model |l^. 

Now we replace the Langevin equation with the following stochastic differential equation 
(SDE), 



dx* = ^-dt 



p 



a* 



dp* = -u{p^*)p*dt* + ^/2D{p^*)dwt, 



(5a) 
(5b) 



Here we have introduced the Wiener process wj. and its difference by 



dwt* = wt*+dt* - wt. = N{t*)dt*. 



(6) 
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The correlations are given by 



{dw;,)KF = 0, (7) 

(rfw"dw^:)RF = dt*5,,5ki. (8) 



The last term of Eq. (l5b 



v/2D(pO*)rfw*. (9) 

is a kind of Stieltjes integral for the stochastic variable. The definition of Stieltjes integral 
in stochastic variables is known to be not unique. Here we consider the three typical cases. 
Note that this problem associated with the discretization scheme is discussed in Refs. 12[ [l^ 
in detail. 

rn 

1. Ito interpretation [17]: 

In this case, the term ([9]) is interpreted as 

v/2^p^ ® dw;, = v/2D(pO*(t*))(w,%,,, - w;). (10) 

Hereafter, we use the symbol ® to indicate the Ito interpretation of the integral mea- 
sure. 



2. Stratonovich-Fisk interpretation 17|: 



In this case (hereafter, we use the symbol o ), the term ([9]) is interpreted as 



= y/D{p^*)dp.y/D{p^*)dt* + V2D(/*) ® dw;,. (11) 

where, in the second line, we used the Ito formula (See Appendix) to convert the 
Stratonovich-Fisk interpretation by the Ito scheme. 



3. Hanggi-Klimontovich 



16 
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20j: 



In this case (hereafter, we use the symbol ^ ), the term as 

^2D{pO*)^dw:, = ^2D{pO*{t* + dt*)){w;,^,,, - w^) 

= 2^/D{pO*){dp^ ^D{p^*))dt* + ^2D{jP*) ® rfw**. (12) 

Here, again, we used the Ito formula to express this interpretation in terms of Ito 
scheme. See Appendix. 
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The equilibrium distribution function described by using these SDEs depends on the 
integral schemes defined above. To see this, we introduce the probability density in the 
phase space defined in the thermal bath rest frame, 

p(x*,p*,f) = (5(3)(x-x*(f))5(3)(p-p*(t*)))Rp. (13) 

The time evolution of p(x*,p*,t*) is given by the Fokker-Planck equation in a unified way 
for the different integral schemes as 

9,.p(x*, p*, n = -Y^ dU^p{^% p\n) + J2 (^(p°*)p'>(x*, p*, n) 

i ^ i 

+ di, (Di--(p°*)a;.D"(p°*)p(x*, p*, f)). (14) 

i 

The values of the parameter a correspond to the different discretization schemes: a = for 
the Hanggi-Klimontovich scheme, a = 1/2 for the Stratonovich-Fisk scheme, and a = 1 for 
the Ito scheme. 

From this equation, we find that the corresponding equilibrium distribution (spatially 
homogeneous and static) is given by 

/ r°(P*) vis) \ 
Peq(x*,p*)ocexp [- j sds^^-a\nD{p\v*))\. (15) 

Since the phase space volume element dV = d^x*d^p* should be an invariant measure un- 
der the Lorentz transformation of the reference frame, we conclude that this equilibrium 
distribution (fT3l) should be a Lorentz scalar function, although the proof is not trivial at 
all 2l|, |22| . As we will see later, we can define the transformation property of the noise by 



using this fact. 



B. Relativistic Brownian motion in a moving frame 

We consider the reference frame which is moving with a constant velocity V with respect 
to the rest frame of the heat bath. We refer to as MF-moving frame. The four-momentum 
dp'^ in this frame is then given by the Lorentz transform of dp*'^ as 



dp^" = A{V)dp*^' 

7(V) /3(\/)n^7(V) 
/3(V)n7(V) 7(V)P||+g 




(16) 
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where 7(V) = — and f3(y) = |V|. The projection operators are defined by P\\ = 

nn^ and Q± = 1 — P\\ with n = V/|V|. In the proceeding calculations, we assume that the 
particle is always on mass-shell = a/p^ + during the whole stochastic process so that 
the stochastic variables, dp^* and dp* are not independent. We then have, 

dp^* = dV^^^T^ = { + (1 - «)^^) ^ 

+Dip'n - ^) } df + vmf^f^ ® d^;.. (17) 

where D'{x) = dD{x)/dx. Substituting this expression into Eq. (fT6|) . we obtain the SDE in 
the MF as 

_ f u{u>^p^MV){p^ - (3{V)py) , Zy(^\r /3(F)m^n ) 



where pv = n^^P and 

'p'-KV)py 



B = j2D{u^p^) , + 7(^) (P° - /3(^K + /3(V^)(np^)) Qx} • (19) 



Here m'^ is the four velocity (7, 7V) , normalized as u'^u^ = 1. The presence of a second rank 
tensor B indicates that the noise is no longer isotropic in the moving frame j23i] . 

The last term is not yet transformed because it contains the noise dw*, which is defined 
only in the RF. We first introduce the stochastic noise which shows the property of the 
Gaussian white noise in the MF as 

(rfWt.)MF = 0, (20) 

(dwl^dwl^JuF = dt*5ij5im, (21) 
.(X)j^p denotes the stochastic average of X in the MF. As was assumed in 



Here the symbo 

previous works 3l-|l3|. we may consider that the correlations of the noise term are Lorentz 
invariant in the rest frame of the particle, and we could write 

{dwt)MF = {dw;,)RF = 0, (22) 
l{p){d^ld^i)MF = l{p*){dwltdwl:)J^F = drSijdku (23) 

where dr is the proper time of the particle, 7(p*) and 7(p) are the Lorentz factors of the 
particle in the RF and in the MF, respectively. They are related through 

7(p) = (A(V)A(-p*))°° = 7(V)7(P*)(1 - (24) 
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where, now 7(V) is the Lorentz factor associated to the Lorentz transformation from the 
RF to the MF. In this case, we conclude that the transformation property of the noise is 
given by 



rfw; = y ^dw, = Vlip)r'iP*)dw, = ^'/\V)^^^^P^d^,. (25) 

Ahhough very reasonable, a proof of the above argument is not obvious. The reason 
is that the Stieltjes integral associated with the noise term is defined on the time span dt, 
so that dw^, is non-local in the time t. Thus the Lorentz transformation entangles with 
the integration scheme in the order of dt. Then the noise term itself is not necessarily 
covariant but can constitute a Lorentz vector only together with the force term. That is, 
the force part and the stochastic part could be mixed. This might be understood well from 
the argument of the kinetic derivation of hydrodynamics. In the rest frame of fluids, the 
velocities of molecules are completely random and all the kinetic energy of molecules are 
replaced by internal thermodynamic quantities. On the other hand, when fluids move, a 
part of velocities of molecules contributes to the collective flow of fluids. In this case, the 
first order correlation calculated in the MF of the noise term in the rest frame of the heat 
bath would not necessarily vanish, 

(dw;)MF ^ 0, (26) 

even if ((iw^,)RF = 0. 

Considering the finiteness of dt, let us write, instead of Eq. ( 125|) . the following more 
general transformation property of the noise, 

jdt* 

dw^, = y —^dwt + Cpdt 



,y^iV),r^J^d^, + Cpdt. (27) 



Here Cpdt term entangles with the force term in the MF, which should be separated from 
the pure stochastic part dwt, satisfying Eqs. ( I20l) and (I2T1) . With this definition and using 
the Ito formula, the Langevin equation is given by 

dp = ^ (1 _ |p _ dt 

+KVhiV)^^ (2 + :^TO^^$(vw) + ^ Cprft + B ® dw„ (28) 
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where 



B = J2(1Q(p!_W>x)b. (29) 



The last term in Eq. ( 128|) should be calculated according to the Ito scheme. Note that the 
Ito formula to convert the SDE from one scheme to the other must be used only after the 
transformation of the noise, Eq. fl27p . 

The corresponding Fokker-Plank equation for the SDE ( l28l) is 



(30) 



where 



, K^>m)7(^)(/-/3(V^)pv) , ^ 



, ,,iV)^ (2 + ^.(^^(^o^V),,). ) n + B C (31) 



Consequently, the equilibrium distribution function is given by the solution of the following 
equation, 

-A^ + ^9^(BB^)^^- peq(x,p) = 0, (32) 
leading to the equilibrium distribution peq satisfying 

{dlPe,)/p., = 2 J2(^B%' - Id', (bB^)'') . (33) 

jkl ^ ^ 

Since peg should be a scalar function, the above expression must coincide with the Lorentz 
transform of the logarithmic derivative of the equilibrium distribution function obtained in 
the RF. A lengthy but straightforward calculation shows that peg can be the Lorentz scalar 
only when 

Cp = 0. (34) 

Thus, for the present model of a Brownian motion with the noise from a given heat bath, 
we conclude that the noise term transforms separately from the force term when we go to 
one reference frame to the other. That is, the transformation of the noise (1251) is consistent 
with the special relativity, and the corresponding SDE in a MF should be given by Eq. ( |T8|l . 
When the background does not satisfy this condition, it is not obvious if we can always 
assume Cp = from the beginning. Our discussion here is only for the thermal background 
is homogeneous and static. 
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C. Generalized fluctuation-dissipation relation 

So far, we have not specified the parameter of the SDE, but when we demand that the 
equihbrium distribution function should be given by the Jiittner distribution in a MF, 

Peq = Const. X e"""*'''/^, (35) 

where T denotes temperature. From this condition, the parameters of the SDE satisfy the 
following relation, 

= ^ (^^^ - «^'(«%)) , (36) 

which is the generalized Einstein's fluctuation-dissipation relation of relativistic Brownian 
motion. Note that this relation depends on a, showing the discretization scheme dependence 
when D is not a constant. 



III. RELATIVISTIC STOCHASTIC ENERGETICS 

Once we have established the SDE for a relativistic particle embedded in a heat bath, 
we can discuss the thermodynamic property of the system defined as the ensemble of these 
relativistic Brownian particles. For this purpose, from now on, we discuss only in the rest 
frame of the heat bath, 

(ix* = ^dt*, (37a) 
dp* = - V0 dt* - u{p°*)p*dt* + ^/2D{pO*) ® dwt* , (37b) 

where we have not assumed the generalized fluctuation-dissipation relation. For the sake of 
the discussion of the first law of thermodynamics, we introduced a scalar potential 0(x*,t*) 
where the explicit time dependence of represents the effect of some external forces. 



A. The First Law of Thermodynamics 

Following Ref. [l|, let us introduce heat as a stochastic variable associated with Brow- 
nian motion. Among the three contributions of the forces in Eg. ( 13 7b I) , the last two terms 
represent the interactions between a system and the heat bath. The first term represents 
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the mechanical force through the scalar potential 0. Thus the work done by the heat bath, 
that is, the heat transfer, d'Q to a Brownian particle is defined as 30 1 



d'Q 



dwt* 
dt* 



o dx*. 



(38) 



Using Eq.f l37b|) we can rewrite Eq.f pSj) as 



d'Q 



dt* 

d{p'* + <P)- 



V0 ) O dx* 

dt* 



dt* 



(39) 



Here, we identify the first term represents the change of the total energy of a Brownian 
particle dE including its potential energy. The quantity {dip/dt*) dt* in the last term is the 
change of the energy contained in the scalar potential due to the change in some external 
parameters contained in 0. That is, it can be interpreted as the mechanical work d'W done 
by the external forces to the system. Thus, rewriting Eq. fl39|) as 



d'Q = dE- d'W, 



(40) 



we can view this as the first law of thermodynamics in the framework of the stochastic 
energetics, the generalization of Ref. |l| for relativistic Brownian motion. Note that this 



relation is satisfied for each event of Brownian motion 



3l| 



B. The Second Law 



In the above, we introduced the potential to identify clearly the role of external work 
in the first law of thermodynamics, and, there, the choice of the parameters of the SDK 
was arbitrary. In the following discussion of other thermodynamic relations, however, we 
will consider only the case where = and impose the generalized fluctuation-dissipation 
theorem to the SDE. 

For deriving the second law of thermodynamics in the stochastic energetics, we introduce 
Shannon's information entropy, 

S=-J dT p{x*,p*,t*)\np{x*,p*,t*), (41) 

where p{x*,p*,t*) is defined in Eq. fll3l) and satisfies the Fokker-Plank equation (fT^ . The 
integration in Eq. (HTl) is taken over all the phase space volume. From these definitions for 
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heat ( 138|) and entropy (HTI) . it is easy to show that the time derivatives of heat and entropy 
are related as 

Since the right-hand side of the above equation is non-negative, we obtain the following 
inequahty, 

which is nothing but the second law of thermodynamics. For the non-relativistic case, see 
Ref. [l|, Q]. An important fact is that, differently from the first law, the second law is 
satisfied only for the average heat, and can be violated for each event, as is the case of the 
fluctuation theorem. 



C. Relativistic Equipartition Theorem 

In the preceding section, we have shown that the first and the second laws of thermody- 
namics can be consistently derived by generalizing the stochastic energetics to relativistic 
Brownian motion. Let us apply this approach to discuss another example, the equipartion 
theorem. 

In non-relativistic systems, the equipartion theorem tells us that the average kinetic 
energy of a particle of an ideal gas in equilibrium is equally distributed by T/2 for each 
degree of freedom. For 3-dimensional mono-atomic gas, we have 

where {0)t denotes the expectation value of O at the instant t, which can be expressed as 

{0), = j rfrp(x,p,t)0(x,p), (45) 

with p(x, p, t) is the solution of the Fokker-Plank equation corresponding to the following 
non-relativistic SDK, 

dx = —dt, (46a) 
m 

dp = -u—dt - V^UhL)dt + V2D ® dwt, (46b) 
m 
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where ?7(x) is a potential and v = D/T. Note that, the above average (H5!) using the solution 
of the Fokker-Planck equation is equivalent to the event average for the corresponding noise. 

In the stochastic energetics, this relation for equipartition of energy is seen to be valid 
in more general situation if the average heat transfer vanishes even if the system is not in 
equilibrium pj]. To see this, we consider the case where D is a constant so that the noise is 
additive. Then, the time derivative of the heat described by this SDE is given by 



dt I niT J \2m 2 
//p^\ _3 



mT\\2m/^ 2 ^ ' ^^^^ 
where the left hand side average is the event average. One can see that the equipartition 
relation is satisfied, when the heat transfer disappears even p is not the equilibrium distribu- 
tion. Thus we may consider the condition of the null heat transfer leads to the equipartition 
relation in more generalized situation. We will apply this idea to relativistic Brownian 
motion. 

The relativistic analogue is not trivial in the sense that there is no simple interpretation 



of the quantity involved 



251]. The expectation value of relativistic kinetic energy, K 



^/p^ + — m does not satisfy a simple relation such as Eq. (jH]). Tolman proposed to read 



Eq. as [261, 



or equivalently. 



where the subscript eq denotes the expectation value at equilibrium, that is, t — > 00. The 
left-hand side of Eq. (HHl) is proportional to the average of p ■ v, which is interpreted to be 
the momentum transfer rate than the kinetic energy itself. 

As an application of the stochastic energetics for relativistic Brownian motion, we can 
derive the relativistic analogue of Eq. P7|) . Taking the noise average of Eq. (l38l) . we obtain 

/ RF \ / RF \ at / j^p 
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To derive this expression, we have used the relation, 



^odp* = dp* ®(^ + U^-) ■ (51) 
From Eq. ( 150|) the null heat transfer rate leads to a relation, 

independent of p. We call the above equation as the generalized equipartition theorem for 
the relativistic Brownian motion. Note that this relation does not have the same form of 
the Tolman relation (149|) . Only when the multiplicative noise is chosen as D oc Eq. ( 152|) 
coincides with 

4\ ='-T, (53) 

as in the non-relativistic case of Eq. fj47|) which is valid even in out of equilibrium. In the 
ultra- relativistic limit, this relation is reduced to 

(E),, = 3T, (54) 

where E = so that it is related to the specific heat of the system. 

The meaning of the equipartition theorem is modified for more general cases. For example, 
suppose the noise has the energy dependence as D oc (p°*)^ with q = const, then Eq. (!52l) 
becomes 

{E^),.=Tx{q + 2){E^-')^^, (55) 
in the ultra-relativistic limit, and 

(|pr+\ =Tmx(g + 3) (|p|%, (56) 

in the non-relativistic Brownian motion (for non-relativistic case, Eq. fl50l) is still the general 
form of Eq. fj47Il when the noise is multiplicative, provided that p^* and D{jp*) are m and 
I?(p^), respectively). These expressions show that the ratio of two moments of kinetic energy 
should be kept constant proportional to the temperature. In this sense we call Eq. ( !52l) as 
the generalized equipartition theorem. 

For q = —2 in Eq. ( l55l) or g = —3 in Eq. (I56l) . the left hand side is finte, but the 
pre-factor in the right hand side vanishes. This means that the average value in the right 
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hand side in each equation diverges but the products 

hm (g + 2) {E'^-')^^ , (57) 
hm (g + 3)(|p|%, (58) 

have finite positive values and the equahties are stih vahd. This is easily checked for equi- 
librium cases, since we can calculate explicitly 

{E'^U = Y^''+'r(g + 3), (59) 
(Ipr^')eq = ^(2mr)(^+^)/2r((g + 5)/2), (60) 

where Z is the normalization factor of the equilibrium distribution function and r(x) is the 
Gamma function. From the recursion relation of the Gamma function, we obtain Eqs. ( I55ll 
and ( l56l) . Of course these relations do not make sense for q < —2 in Eq.(l55l) or for q < —3 
in Eq.f l56|) . since the left had side (so the right hand side, too) diverges. 

Note that the condition, Eq.( l55|) or Eq.( l56|l . itself does not determine the distribution 
function. Of course, if we assume that there exists a distribution function, f{E;t*), inde- 
pendent of q, satisfying Eq. (l3^ for any g > 2, we can demonstrate that such a distribution 
function should be proportional to exp{—E/T), and similary for non-relativistic case. How- 
ever, for a finite t*, the distribution function still refiects the initial condition and depends 
on q. 

Thus, as stated before, the null heat transfer condition does not mean that the system is 
in equilibrium, but inversely, the null heat transfer condition is naturally satisfied in thermal 
equilibrium for any choice of D{p^*). For more general noise coefficient, this can be seen 
directly from Eq. ( I50l) by re-expressing it as 

f - /"^ w ^ (^) ■ '^''-<^---''-'f 

(61) 

The integrand vanishes if p is given by the Jixttner distribution function. Therefore, in 
thermal equilibrium, both of the generalized equipartition theorem f l52|) and the Tolman 
relation ( l49ll are satisfied at the same time. 

To see how the null heat transfer condition is attained, we consider the case of a pion 
m = 139 MeV in a heat bath of T = 150 MeV as an example for a possible application to 
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FIG. 1: The two ratios as a function 
of time t* with the initial condition 
of 50 MeV. 

RIC physics. In Figs. [H and El we plot the ratio, 
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FIG. 2: The two ratios as a function 
of time t* with the initial condition 
of 300 MeV. 



Rqe {i* 



T 



g(p°*)p* 



(62) 



t* 



by solving Eq. (]37b|) for one dimensional case with = 0. In this example, we compare 
the two different situations; one is that the pion is initially has a lower energy than the 
temperature of the heat bath, p(t* = 0) = 50 MeV (Fig. [1]), and the other the pion has 
initially a higher energy, p{t* = 0) = 300 MeV (Fig. |2]). When the null heat transfer 
{d'Q/dt*)^p = is attained, we should have Rqe (t*) = 1. For the sake of comparison, we 
show also the ratio 



p 



(63) 



which corresponds to the Tolman relation for one dimensional case. In thermal equilibrium, 
this value should stay at -RtoI (t*) = ^- In these simulations, we consider the case of g = — 1, 
that is, D{p^*) = K/p^* with = 0.1 (MeV)^. As seen from these figures, the null heat 
transfer condition is almost satisfied before t* ~ 0.5 fm/c, whereas the Tolman relation only 
converges for t* > 1 fm/c. We checked some different parameter sets and found that such 
a behavior always appears. 
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IV. SUMMARY AND CONCLUDING REMARKS 

In this paper, we extended the stochastic energetics for relativistic Brownian motion. 
As a model for the relativistic particle embedded in a heat bath, we first established a 
relativistic SDE which produces the invariant Jiittner equilibrium distribution under an 
arbitrary Lorentz transformation. Using this model, we discussed the thermodynamic laws 
and the equipartition theorem, by applying the stochastic energetics. We showed how the 
first and second laws of thermodynamics are derived from relativistic Brownian motion in 
this context. We obtained the explicit form of the heat transfer rate between the relativistic 
Brownian particle and the heat bath. Of course, our results recovers the corresponding 
results of Ref. 1[ in the non- relativistic limit. As a result, we showed that the concept of the 
stochastic energetics is applicable to relativistic Brownian motion. 

We further showed that the condition of the null heat transfer leads to the generalized 
equipartition theorem. Except for the particular choice of the parameters of the SDE, the 
generalized equipartition theorem does not coincide with Tolman's relativistic equipartition 
theorem if not in thermal equilibrium. The null heat transfer and thermal equilibrium are 
not equivalent, but the former includes the latter. As a matter of fact, we found that 
the generalized equipartion relation is attained before the system reaches the true thermal 
equilibrium in a few examples. It will be interesting to investigate more in detail the meaning 
of null heat transfer condition with respect to possible transient thermodynamic properties. 

The stochastic energetics is considered as a promising approach for the study of thermo- 
dynamics of mesoscopic systems jll, Q]- Thus, the present study will be a starting point of 
the application of the stochastic energetics for the physics of relativistic heavy ion collisions. 
There, the finite size of the system, as well as the short reaction time are not negligible so 
that the deviation from the local thermal equilibrium must be clarified. In particular, the 
approach via the stochastic energetics allows introducing the energy conservation between 
microscopic and macroscopic degrees of freedom. This may serve for a consistent descrip- 
tion of the formation of initial condition or the freezeout processes of final states in the 
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Appendix A: Ito formula 

Let us consider an arbitrary function /(x) and the evolution of x is given by the SDE, 

c?x = Adt + B ® dw. (Al) 

Then, the variation of /(x) is 

d/(x) = {Y, A%f{^) + i J^lBB^f^ didjf{^)}dt + B*^9./(x) ® dw^. (A2) 

i ij ij 

This is called the Ito formula 17|. By using the Ito formula, we obtain. 



[B o dwf = [B ® dwf + ^Y1 ^^''9jB'''dt, (A3) 

jk 

[B * dw]' = [B® dwf + J2 B^''djB'''dt. (A4) 

jk 

Thus we can conclude as follows. When we have the Stratonovich-Fisk SDE, 

rfx = Adt + B o dw, (A5) 

this is equivalent to the Ito SDE, 

(ix* = {A^ + ^Y1 ^^'"djB'''}dt + B®dw. (A6) 

jk 

When we have the Hanggi-Klimontovich SDE, 

dx = Adt + B dw, (A7) 

this is equivalent to the Ito SDE 

dx.' = {A' + Y B^''djB'''}dt + B®dw. (A8) 

jk 
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